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ABSTRACT 



Several inflationary models predict the possibility that the primordial perturbations of 
the density field may contain a degree of non-Gaussianity which would influence the subse- 
quent evolution of cosmic structures at large scales. In order to study their impact, we use a 
set of three cosmological DM-only simulations starting from initial conditions with different 
levels of non-Gaussianity: /nl= 0,±100. More specifically, we focus on the distribution of 
galaxy clusters at different redshifts and, using suitable scaling relations, we determine their 
X-ray and SZ signals. Our analysis allows us to estimate the differences in the logA^-logS' and 
logA^-logF due to the different initial conditions and to predict the cluster counts at different 
redshifts expected for future surveys (eROSITA and SPT). We also use a second set of sim- 
ulations assuming a different cosmological scenario to estimate how the dependence on /nl 
is degenerate with respect to other parameters. Our results indicate that the effects introduced 
by a realistic amount of primordial non-Gaussianity are small when compared to the ones 
connected with current uncertainties in cosmological parameters, particularly with Og. How- 
ever, if future surveys will be associated with optical follow-up campaigns to determine the 
cluster redshift, an analysis of the samples at z > 1 can provide significant constraints on /nl- 
In particular we predict that the SPT cluster survey will be able to detect 1000 clusters at 
z > 1 for the Gaussian case, with a difference of 15-20 per cent associated to /nl= ±100. 

Key words: cosmology: theory - (cosmology:) large-scale structure of Universe - galax- 
ies: clusters: general - X-rays: galaxies: clusters - methods: A^-body simulations - methods: 
statistical. 



1 INTRODUCTION 

The universally accepted scenario for the formation of cosmic 
structures in the universe is based on the mechanism of gravita- 
tional instability, which assumes that the density fluctuations gen- 
erated at some early epoch grow by accreting mass from the sur- 
rounding regions through gravitational processes. The origin of 
these cosmological seeds is generally related to the final phases 
of the inflationary expansion for which a l arge variety of theoreti- 
cal models exist in the literature (see , e.g., lKi nnevll2008l : iLangloi j 
l2008l : lBaumann & Peirisl2008l : lBaum ann 2009, for recent reviews). 
These models originate perturbations having different statistical 



properties, usually investigated in terms of probability distribu- 
tion function (PDF) and correlation functions/power polyspectra. 
In particular, the most standard slow-rolling models, where a single 
field is responsible for the inflationary accelerated expansion, pro- 
duce fluctuations having almost uncorrelated phases. For this rea- 
son in cosmological studies it is usual to assume that the primordial 
perturbations are Gaussianly distributed, which leads to the further 
simplification that their complete description is possible using the 
power spectrum only. However, even the simplest inflationary mod- 
els allow for small departures from Gaussianity, which can become 
more significant in non-standard models, like the scenarios based 
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on the curvaton, the inhomogeneous reh eating and the Dirac-Born- 
Infeld inflation (see lSartolo et al.l2004l and references therein). As 
a consequence, the observational determination of the amount of 
non-Gaussianity present in the primordial fluctuations is now con- 
sidered not only a general probe for the inflationary concept, but 
also a powerful discriminatory test between its various theoretical 
models. 

Given the infinite variety of possible non-Gaussian models, 
it is necessary to introduce a simple way to quantify the level of 
primordial non-Gaussianity. In the recent years it has become stan- 
dard practice to adopt the dimensionless non-lineari t y parameter 
/nl ( see Salopek&Bond 199G; Gangui et al.lll994 IVerde et al.1 
120001 : iKom atsu & Spergel 2001), measuring the importance of the 
quadratic term in a sort of Taylor expansion of the Bardeen's gauge- 
invariant potential namely 



pendent probe for /nl, in particular they showed that the most 
evident non-Gaussian signatures are present in the mass func- 
tion and clustering (bias and bispectrum) of dark matt er mass 



* = *L+/NL(*t-(*L)) 



(1) 



here <&l represents a Gaussian random field. In particular, in this 
paper we follow the so-called large-scale structure (LSS) conven- 
tion, where <I>, that on scales smaller than the Hubble radius corre- 
sponds to the usual Newtonian peculiar potential (but with changed 
sign), is linearly extrapolated to the present epoch. 

For many years, the study of the statistical properties of the 
cosmic microwave background (CMB) has been considered the 
most efficient way to measure /ml- In fact its temperature fluctu- 
ations are directly related to the density perturbations in a regime 
in which the non-linearities originated by the subsequent process 
of gravitational instability are not modifying their primordial char- 
acteristics, including the PDF. Different statistical estimators have 
been applied to the most recent data, keeping improving constraints 
on /nl- So far the more stringent results come from the analyses 
of the 5-years WMAP dataset: assuming a local shape for non- 
Gaussianity (as we will do in this paper), iKomatsu et alj ^20091) 
foun d that /nl varies between -12 and 145, while Smith et aU 
12009) found —5 < /nl < 104Il(s ee also the positive dete ction of 
non-Gaussian features reported bv lYadav & Wandei3l2008h . 

More recently it became clear that the LSS represent an al- 
ternative tool, potentially as valid as the CMB to constrain /nl- In 
fact, deviations from an initial Gaussianity induce a different tim- 
ing for the whole process of structure formation, providing an in- 
teresting framework to look for specific non-Gaussian imprints. In 
general, as already evident from the results of the firs t generation of 
non-Gaussian Af-body simulations in the early '90s jMessina et al.l 
Il990l : lMoscardiniet"ai]|l99ll : [Weinberg & Coldl 19921) . if the pri- 
mordial density field is positively (negatively) skewed, the for- 
mation is favored (disfavoured) and structures of a given mass 
form at earlier (later) epochs. However, there is an important 
difficulty in using LSS for constraining /nl: the late non-linear 
evolution introduces additional non-Gaussian features that need 
to be disentangled from the primordial ones. For this reason it 
has been necessary to have resort to suitable high-resolution A'- 
body simulations to follow the growth of the LSS also in the 
full non-linear regime, and to calibrate the expected signatures as 
a function of the primordial non-Gaussianity. This has been ex- 
tensively done in the recent years by differen t groups ( see, e.g 




Rang et alj|2007l: iGrossietd] l200l [Palal et al. 2008; Viel et al 



Desiacaues et al.1 l2009l: iPillepich et Zl [20081 : IGrossi et al 



The results allowed to assess the power of LSS as inde- 



' We multiplied by a factor of 1.3 the original results to convert them to the 
/nl LSS-convention adopted here. 



haloes (see also the analytical predictions made by Matarrese e t al 
2000: Lo Verde et al. 2008; Afshordi & ToUey 2008; Carbone et al 
2008; McDonald 2008; Maggiore & Riotto 2009; Lam & Shett 



,2009; Valageas 2009; Verde & Matarrese 2009). The first attempts 
of an application to real data gave very encouraging constraints: 
[siosar et al. ( 2008), combining the bias measurements for two sam- 
ples o f luminous red galaxies a nd quasars, found /nl = 37^^^, 
while lAfshordi & TollevI ( |2008|) . studying the integrated Sachs- 
Wolfe effect (ISW) in the NVSS survey, derived /nl = 272 ± 127; 
error bars are at 2-o level. Notice t hat in both cases we report the 
values as revised by iGrossi et alj ( l2009t) to include a correction 
mimicking the ellipsoidal collapse. 

Being at the top of the hierarchy of structure formation, galaxy 
clusters are in principle ideal probes for constraining /nl- Indeed, 
the statistics of fare events, either galaxy clusters or deep voids 
jViel et al.|[2009l) are very sensitive to primordial non-Gaussianity. 
However, until now the use of clusters as non-Gaussian probes has 
been hampered by two practical problems: first, the observational 
estimate of their mass is affected by large uncertainties, whatever 
are the method and the observational band adopted; second, it is 
difficult to build up samples that span a large range of redshift and 
are also statistically complete. The first problem can be in some 
way overcome by using the scaling relations (expected from theo- 
retical arguments and confirmed by observations) existing between 
mass and different observables (see below). The second problem 
will be solved in the next years, thanks to the efforts of a set of 
(in progress or planned) surveys, which promise to provide large 
samples of galaxy clusters covering a volume comparable to the 
horizon size: see, e.g., DBS, PanSTARRS, BOSS, LSST, ADEPT, 
EUCLID. 

In this paper we will focus on estimating the signals pro- 
duced by galaxy clusters jn the X-ray band and through the 
Sunyaev-Zel'dovich (SZ, fSunvaev & Zel'doviclj|l97j|. Il98d) ef- 
fect in different non-Gaussian scenarios. We will take advantage 
of the fact that, in both cases, the corresponding observables, i.e. 
the X-ray luminosity and the Compton j-parameter, are related 
to the mass by well calibrated scaling relations. This will allow 
us to introduce selection criteria mimicking the characteristics of 
specific surveys. In particular we will c onsider the properti es of 
the wide surveys plarmed w ith eROSITA jPredehl et al]|2007l) and 
SPT (' Carlstrom et^l2009b . as examples of future X-ray and SZ 
projects, respectively. The main goal of this work is to figure out 
what are the observational evidences of the presence of some level 
of primordial non-Gaussianity, discussing the possibility of con- 
straining /nl with these future datasets. Notice that the same sam- 
ples h a ve bee n considere d for a s imilar work in ISefusatti et al.l 
bood) : lOguril feoOft) : Fedeli et alj ( |2009) . In this paper we ad- 
dress this problem using numerical rather than analytical tools. The 
advantage is twofold. First of all, A'-body simulations permit to 
fully account for non-linear evolution which, instead, is usually 
accounted for by analytical models in an approximate way only. 
Second, numerical experiments allow us to extract realistic mock 
cluster catalogs that can be easily used for modeling the observa- 
tional selection, which is more difficult to account for in a Fisher 
matrix-like approach. 

The plan of this work is as follows. In Section |2] we present 
the numerical simulations of non-Gaussian models on which the 
following analysis is based; we also describe the method applied 
to construct the light cones. Section[3]introduces our model for the 
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X-ray emission from galaxy clusters and reports the corresponding 
results in terms of number counts, paying attention to the expected 
results for the eROSITA wide survey. Section |4] is devoted to the 
model for the SZ signal and to the corresponding results, given in 
terms of source counts and statistical properties of the maps; the 
specific case of the SPT wide survey is treated. Section[5]discusses 
the possibility of using the differential redshift distribution to con- 
strain the primordial /nl- Finally, in Section[6]we draw our conclu- 
sions. 



2 MODELS AND METHOD 
2.1 The simulation sets 

In order to study the possibility of detecting the signatures of pri- 
mordial non-Gaussianity in the LSS of the Universe we must take 
into account the complete process of structure formation. For this 
we make use of the outputs of two different sets of DM-only simu- 
lations, focusing on the distribution of the DM haloes associated to 
the galaxy clusters, as identified in the dif ferent snapshots. 

The first simulation set (presented in lGrossi et a l. 2009) con- 
sists of three DM-only simulations performed with different lev- 
els of primordial non-Gaussianity, that, expressed in terms of the 
dimensionless non-linearity parameter /nl, are /nl=0, ±100; the 
case /nl=0 corresponds to the standard Gaussian case. The ini- 
tial conditions were set by assuming a flat ACDM model dom- 
inated by a cosmological constant with para meters chosen to be 
consistent with the WMAP three-year results jSpergel et ai]|200l 
WMAP -3 hereafter): namely, the density contributions of cosmo- 
logical constant and matter correspond to £1^=0.76 and fl,n=0.24, 
respectively, while the normalization of the power spectrum of den- 
sity fluctuations is fixed as 08=0.8, being Og the r.m.s. matter fluc- 
tuation into a sphere of radius 8/i^' Mpc. The three simulations 
started from the same random generation of initial conditions with 
the only difference consisting in the different value o f /nl- All the 
runs, carried out wit h the N-hoAy code GADGET-2 ISpringel et al.l 
l200ll ; ISDring"eill2005h . followed the evolution of 960^ DM particles 
inside a cubic volume of 1200/i^' Mpc per side, with each parti- 
cle having a mass of m = 1.4 x lO^'/i^'M©; here h represents the 
Hubble parameter defined as, h = HqI(\QQ km s^' Mpc^')=0.7. 
The gravitational force has been computed using a Plummer equiv- 
alent softening length £ = 25/;^' kpc. 

These simulations produced 14 outputs in the redshift range 

^ z ^ 4. For each snapshot we produced a catalogue of DM 
haloes identified using a.friends-of-friends (FOF) algorithm, adopt- 
ing a linking length of 0.2 times the mean interparticle distance: 
with this choice the size of the haloes roughly corresponds to their 
virial mass, Mvii-. Since in this work we are interested in the X-ray 
and SZ signals produced by galaxy groups and clusters, we only 
considered DM haloes having Mvir > lO'^'/i^'My,. For each halo 
we kept the information on the mass and the position of its center 
inside the cosmological volume. 

The second set of simulations is the one described in 
iGrossi et al. I j2007t) . It consists on 7 runs covering a wider range 
of non-Gaussianity: /nl=0, ±100, ±500, ±1000. The main dif- 
ferences with respect to the first set are that the box size is only 
500/;^' Mpc with 800^ particles and that the c osmological model 
adopted is close to a WMAP first-year jSpergel et al. 2003, WMAP- 

1 hereafter) cosmology, with the following cosmological parame- 
ters: O.\=0J, 0.^=0.3, ©8=0.9. For this set we have 21 outputs in 
the range ^ z ^ 4 and a halo catalogue for each snapshot ob- 
tained in the same way as described before and with the same mass 



limit. Notice that in this paper this last set of simulations will be 
mainly used to test the impact of different cosmological param- 
eters compared to non-Gaussianity: for this reason we will dis- 
cuss only the results of the most extreme and the Gaussian model 
(/nl= 0, ±1000). This simulation set will also be used in order to 
check the effect of the finite box size on our results. 

2.2 Constructing the light-cones 

As said in Section [T] in order to study the impact of primordial 
non-Gaussianity on the LSS we adopt an observationally-oriented 
approach. To this purpose we use the halo catalogues described in 
Section [ZT] to produce mock light-cones by stacking several sim- 
ulation volumes. In particular we want to cover the redshift range 
^ z ^ 4, which corresponds to a comoving distance of 5249 h^^ 
Mpc (5019 /i^' Mpc with the cosmology adopted for the second 
set). This length requires to stack 5(11 for the second set) times the 
simulation box. However, in order to obtain a better redshift sam- 
pling, we divide the simulation volume into slices along the line 
of sight. The number of slices varies from cube to cube and their 
comoving distance intervals are created in order to allow us to use 
all of the 14 (21 for the second set) snapshots. More precisely, for 
any given distance from the observer we compute the correspond- 
ing time elapsed from the big bang and we choose the snapshot that 
better approximates this value. 

In order to avoid the repetition of the same structures along 
the line of sight, for every stacked simulation volume we perform 
a randomization of the halo spatial coordinates: for every cube we 
randomly choose the axis to put along the line of sight, we assign 
a 50 per cent probability to reflect each axis and, since our simu- 
lations assume periodic boundary conditions, we proceed to a ran- 
dom recentering of the coordinates. In order to preserve the whole 
information on the structures inside the simulations' volume, the 
slices belonging to the same cube undergo the same randomization 
process. With th i s met hod, which is similar to the one adopted by 
iRoncarelli et al. we obtain 18 different slices belonging to 

5 independently randomized cubes (3 1 and 1 1 for the second sim- 
ulation set, respectively). This process is repeated with the same 
initial random seed for all the simulations of the set. Each light- 
cone produced in this way spans an angle of 13.1 deg (5.71 deg for 
the second set) per side, determined by the length of the box at the 
maximum redshift, z = 4. 

By varying the initial random seeds we can obtain different 
light-cone realizations that we can use to assess the statistical ro- 
bustness of our results. For each non-Gaussian model, we created 
20 (100) different light-cone realizations, thus covering a total area 
of 3432 deg^ (3260 deg^). However, it is important to note that 
this area cannot be considered as completely independent, being 
produced starting from the same finite volume of the simulation: 
as a reference, at z ~ 4 the same simulation volume is completely 
stacked in all realisations, while at z — 0.5 we are able to produce 
about 15 independent volumes crossing the light-cones. 

Using the whole set of light-cone realizations we compute the 
mass functions, expressed in terms of number of objects per solid 
angle, for 4 redshift bins. We show in Fig.[T]the results for the three 
models of the first set (/nl=0,±100) and for the /nl=0,±1000 
models of the second set. 

Looking at the ^ z ^ 0.5 bin, the three mass functions of 
the first set are basically indistinguishable. The non-Gaussian de- 
viations become more evident in the tail at higher masses and at 
higher redshifts: in the second bin (0.5 ^ z ^ 1) the /nl= ±100 
models introduce a ±5 per cent difference in the high mass (Myh- > 



4 M. Roncarelli et al. 



1000.0 
100.0 
10.0 
1 .0 
0.1 






-0 
-0 







1 000.0 




0.0<z<0.5 : 


100.0 






10.0 






1 .0 






0.1 




0.4 
0.2 
0.0 
-0.2 
-0.4 





1 3.0 

1 000.0 
100.0 
10.0 





-0 
-0 



0.4 
" .2 
.0 
.2 
.4 



13.0 













0.5<z<1.0 - 


















■ - 







1 3.5 



1 4.0 



14.5 

















1.0<z<2.0 : 















15.0 13.0 

1000.0 
100.0 
10.0 
1 .0 






.1 





4 





2 











2 





4 



1 3.5 



14.0 



1 4.5 



1 5.0 

















2.0<z<4.0 : 












— ' — \ — ' — 


~x "\ "x 
\ ~x 

^\ ■\ \ 


fNL= 100 ^ 

fNL= ; 

fNL= -100 , 






' a/^ " ' 





























4.5 



1 5.0 



13.5 14.0 14.5 15.0 13.0 13.5 14.0 

Log(M,[yh"^Mo) Log(M„;yh"^MQ) 

Figure 1. The logW-logM for the two simulation sets in different redshift bins. In the upper panels the solid lines refer to the first simulation set (average of 
20 light-cone realizations) for the three different levels of non-Gaussianity: /nl=0,±100. The dot-dashed lines refer to the second simulation set (average of 
100 light cones, WMAP-1 cosmology) for models with /nl=0, ±1000. In the lower panels we show, for the first set of simulations only, the difference AA? /N, 
computed with respect to the Gaussian (/nl=0) simulation. 



10 /j- Mq) cluster counts, which increases to about 10 per cent 
at z > 1. This is in agreement with the fact the modification of the 
distribution of the primordial density fluctuations primarily affects 
the form ation of the biggest haloes at early epochs, as already dis- 
cussed in lGrossi et alj ( l2009h . 

When comparing the two simulation sets, it is clear that the 
differences between the various non-Gaussian scenarios are very 
small compared to the ones resulting from the change of the cos- 
mological model. For example, even at z > 1 the /nl= 100 model 
of the first set adds only about 10 per cent counts to the Gaussian 
case (Mvir > IO^'^/j^'M©), while the WMAP-1 cosmology scenario 
predicts more than three times as much objects. This fact highlights 
that the uncertainties with cosmological parameters, particularly 
with Og, are critical when addressing the problem of detecting pri- 
mordial non-Gaussianities in the L SS. More in the detail, an alytical 
models of the mass functions (e.g. ISheth & Tormenir2002h predict 
that a difference of ±0.01 in Og produces a change of about 10 per 
cent in the total integrated counts (Mvii- > IO^^/i-'Mq), raising to 
about 15 per cent at z > 1, where the /nl parameter is expected to 
produce the most significant effects. This aspect will be considered 
in the following analyses. 

It is worth to note also that for masses higher than ~ 
W^'^-^h-^MQ the mass functions of the second simulation set 
steepen and approach the first-simulation ones: this is an artificial 
effect produced by the smaller box size of the second simulation 
(500 instead of 1200/;"' Mpc). In fact this causes a loss of power 



for perturbations at large scales, corresponding to a lower abun- 
dance of objects with high masses, and gives an indication of the 
range of validity of our simulation sets. 



3 MODELING X-RAY COUNTS 

Since our simulations consider only DM particles, in order to com- 
pare our results with present and future cluster surveys we need 
to define a model to associate the baryonic component to each 
DM halo. In particular we will focus here on modeling the X-ray 
emission (and the SZ signal, described in Section|4j, exploiting ob- 
served and predicted scaling relations. As said in Section [ZT] our 
method provides the virial mass Myir of our clusters, while often 
scaling relations are published using other mass definitions (like 
M2oo> ^500 Md so on): therefore, in order to convert Myir to the re- 
quired value, we assume tha t our DM haloes follow a NFW density 
profile jNavarro et al.l 19971) with a concentration param eter c given 
by the c[M,z) relation proposed bv lDolag et al.l l l2004l) . Notice that 
here we explicitly neglect the presence of any diffuse emission from 
the IGM associated to the warm-hot intergalactic medium (WHIM) 
and concentrate only on signals coming from galaxy clusters and 
groups. Even if the WHIM component is expected to contribute 
significantly to the total LSS signal, its presence i s not expected to 
signifi cantly affect the clusters count rate (see, e.g.. lRoncarelli et aH 
l006ah . 
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In order to associate an X-ray luminosity to a DM halo we 
assume the phenom enological mass-luminosity relation found by 
IStanek et alj ( l2006h . namely 

h-^^ loXrgs-i =^'^'°^(^>' ( 101?-1M J ' 

where Lx is the luminosity in the [0.1-2.4 keV] band and M200 
is the mass of the cluster inside a radius enclosing 200 times the 
critical density Pc(z) of the Universe at the redshift of the cluster. 
The normalization Ljj g corresponds to the luminosity of an object 
with M200 = lO'^/j^'Mrr, at z = 0. The term E{z) represents the 
redshift evolution of the Hubble parameter, 

E(z)=a^{\+zf +O.A (3) 

while the value of s is assumed to correspond to the self-similar 
evolution case, s = 7/3. The best-fit p arameters ^15 n and p depend 
on the assumed cosmological model. IStanek et al.l ( |2006[) publish 
their results assuming both i2m=0.24 (InLij q = 1 ■ 19, p = 1 .46) and 
i2rn=0.30 (lnLi5 = 1.34, p = 1.59), therefore we can take those 
values as a reference for our first and second simulation sets, re- 
spectively. We also take into account a scatte r in the Lx—M relation 
that w e fix to 17 per cent, as measured bv iReiprich & BohringeJ 
( l2002h . 

In this work we will determine the observed clusters fluxes 
in the [0.5-2 keV] and [0.5-5 keV] bands to compare our results 
with the abundanc es derived from ROSAT X-ray clusters survey 
jRosati et alj2002l) and to predict the expected counts for the forth- 
coming eROSITA survey jPredehl et aklhoOTl) . respectively. In or- 
der to calculate the band corrections, we need to assume the ICM 
temperature that determines the spectral distribution of the emitted 
radiation. For this purpose we consider our haloes as isothe rmal 
and we use the A/200 ~ ^200 relation which lArnaud et al] i2005h ob- 
tained from a sample of ten nearby clusters observed with XMM- 
Newton, adding a self-similar redshift evolution. Notice that their 
sample covers the temperature range 2-9 keV, so we are forced to 
extend this relation to smaller objects, where a steepening o f the 
^^200 ~ 7*200 relation is expected. In fact, lArnaud et al.l ( |2005|) ob- 
tain a slope higher of ~ 0.2 when restricting their sample only the 
hottest (T > 3.5 keV) clusters. Even if neglecting this effect may 
lead to an underestimate of the cluster temperatures for the small 
objects, we checked that it does not have a strong impact on our 
final results. For example, if we assume a further steepening of 0.2 
in colder {T <2 keV) clusters, this leads to a change of about 1 per 
cent in the faint-end of the logW-logS, and the relative differences 
between the non-Gaussian scenarios remain unchanged. 

It is known that the ICM cools mainly via bremsstrahlung 
emission which is the main physical process responsible of the X- 
ray emission of galaxy clusters. Therefore we can model the emis- 
sion of the gas assuming a free-free spe ctrum with a G aunt fac- 
tor g(£/fcB) 7200 = (^/fcfi)""'^ (see, e.g.. lBorganietalj |T999. for 
more details). With this simplifying assumption we are neglecting 
the presence of other known emission processes like line-emission 
from metals, which can give a non-negligible contribution espe- 
cially to low-temperature (T < 2 keV) clusters; however, for the 
reasons explained above, neglecting this process has no significant 
impact on our final results. 

Once obtained the band correction /band^ we calculate the clus- 
ter flux in a given band as 

where di{z) is the luminosity distance of the cluster. 
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Figure 2. The number of haloes a.s a function of the flux limit in the [0.5-2 
keV] band for the two simulation sets. The solid lines refer to the three mod- 
els of the first set (average of 20 light-cone realizations) with /nl=0,±100. 
The dot-dashed line refer to the second simulation set (average of 100 light- 
cones, WMAP-1 cosmology) for models with /nl=0, ±1000. T he diamonds 
with errorbars correspond to the number counts derived by iRosati et alj 
1 2002.) from ROSAT X-ray clusters survey. 

Fig.|2]shows the number counts as a function of the X-ray flux 
limit in the [0.5-2 keV] band for the three models of the first set 
and for the /nl=0,±1000 of the second simulation set. Our results 
are compared with the counts obtained from the ROSAT survey, in 
the same band. All of the first three models show a good agreement 
with the data, having error bars much wider than the difference 
associated to the change of the /nl parameter. On the contrary, 
the counts obtained with the WMAP-1 cosmology are significantly 
higher (even when considering the standard Gaussian case): the 
only way to fit the data with these cosmological parameters would 
be to assume a strong negative evolution of the luminosity with red- 
shift, e.g. adding an extra factor (1 +2)!^ t o equation[2l with |3 ~ —2, 
which is absolutely unrealistic (see, e.g.. lEttori et al.ll2004l) . These 
results show that the dependence of the cluster counts on different 
cosmological parameters, and particularly on Og, is much higher 
with respect to the one on the value of /nl: this highlights the dif- 
ficulty of constraining the level of the primordial non-Gaussianity 
with the current uncertainties in the cosmological parameters. 

Fig. [3] shows the logA'-log5 in the [0.5-5 keV] band in dif- 
ferent redshift intervals. Again, the number counts differences be- 
tween the different non-Gaussian scenarios are about one order of 
magnitude smaller than the corresponding differences obtained by 
increasing the value of Og from 0.8 to 0.9. Anyway, it is interesting 
to note that in the redshift interval 0.5 < z< 1, the number counts 
expected assuming the eROSITA detection limit (3.3 x 10^'^ erg 
s^' cm^^) differ by about 10 per cent when assuming /nl= ±100 
with respect to the Gaussian model. These relative differences grow 
to about 20 per cent in the 1 < z < 2 interval and increase for higher 
fluxes: this result is in agreement with the expected evolutionary 
scenario emerging from Fig. [T] with higher mass (and then more 
luminous) haloes forming earlier in models with positive /nl- 

Table \T\ shows the number of galaxy clusters with M200 > 
10^"^ h-^M.,, that are expected to be detected by the eROSITA wide 
survey under our assumptions. For our predictions we assume an 
effective area of 20000 deg^ and a detection limit of 3.3 x 10^ erg 
s^' cm^^ in the [0.5-5 keV] band. For the Gaussian model we pre- 
dict ~ 60000 detections, one per cent of them at redshift larger than 
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unity. The changes due to the presence of a moderate level of pri- 
mordial non-Gaussianity (/nl= ±100) are small, but always larger 
than the expected poissonian error: considering the total counts, the 
variation is about 4 per cent, while for the redshift bin 1 < z < 2 the 
abundances change by 15 per cent. However, once again, the de- 
pendence on /nl is much smaller when compared to the one on the 
power spectrum normalization Cg: at z > 1, the differences in the 
cluster counts associated to the /nl= ±100 models are comparable 
to an uncertainty of ±0.01 in ag (see the discussion in Sect. 12.2b . 
For what concerns the total integrated counts, the number of objects 
detectable by the eROSITA survey predicted by the Gaussian simu- 
lation of the second set (Og = 0.9) grows by a factor of 2: ~ 120000 
objects, with a similar redshift distribution. This confirms the ne- 
cessity of having alternative derivations of the main cosmological 
parameters to allow to exploit the power of cluster coun ts to con- 
strain the level of primordial non-Gaussianity (see also iFedeli et al.l 
I2OO9I) . 



4 THE THERMAL SZ EFFECT 

Another important observable quantity to study galaxy clusters is 
the thermal SZ (tSZ) effect, namely the inverse-Compton scattering 
of the CMB photons caus ed by the electrons present in the hot intr- 
acluster plasma (see, e.g.. lBirkinshaw|[l999l : ICarlstrom et al |2002|. 
for detailed reviews). This effect causes a distortion in the CMB 
blackbody spectrum, whose intensity in a given direction can be 
expressed in terms of the Compton y-parameter defined as 



Table 1. The number of detected haloes (with M200 > 10 A" M,; ) pre- 
dicted for the eROSITA wide survey in different redshift bins. The quoted 
errors ai'e assumed to be poissonian. 



Number counts 
/nl < z < 0.5 0.5 < z < 1 1< z < 2 



Total 



-100 41590±204 15934±126 478± 22 58007±241 

42318±206 16715±129 571± 24 59610±244 

100 43443±208 18015±134 641± 25 62105±249 



f 



(5) 



here is the Boltzmann constant, 07- is the Thompson cross sec- 
tion, nie is the electron mass, c is the light speed, while ng and 
Te represent the electron number density and temperature, respec- 
tively. The distortion of the signal results in a difference AT in 
the measured temperature which depends on the observational fre- 
quency. In the Rayleigh- Jeans (RJ) limit this is given by 



-2v , 



(6) 



AT 

TCMB 

where rcMB=2.726 K is the CMB temperature l lMatheretal.ll994t) . 
thus producing a temperature decrement, which can be as high as 
AT ~ 10^^ K for the central regions of the most massive clusters. 
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Table 2. The number of detected haloes (with M200 > IO'^/i^'Mq) pre- 
dicted for the SPT wide survey in different redshift bins. The quoted errors 
are assumed to be poissonian. 



/nl 



Number counts 
0<z<0.5 0.5<z<l l<z<2 



Total 



-100 7796±88 5961± 75 960± 31 14456±120 

7958±89 6041± 78 1124± 34 15137±123 

100 8144±90 6357± 80 1268± 36 15784±126 



Following an observationally oriented approach, the interest- 
ing quantity to be evaluated for a given halo is the integrated Comp- 
ton y-parameter Y, defined as 



; [ ydQ. -- 
Jn 



1 



(7) 



where Q. is the solid angle subtended by the cluster and V is its 
physical volume. This adimensional quantity depends on the an- 
gular diameter distance d^iz), and constitutes an equivalent to the 
flux in the X-rays. Therefore it is useful to introduce the intrinsic 
Compton y-parameter defined as 



Y'"' = YdA(zr 



(8) 



which is roughly proportional to the mass and to the temperature of 
the object. Using hydrodynamical simulations it has been possible 
to calibrate scaling re lations between the cluster mass and its SZ 
observables (see, e.g.. lDiaferio et alj|200^ ; Shaw et al ] |2008l) : here 
we adopt the M — F"" relation, found by Nagai ( 2006h . that can be 
expressed in the form 



-•200 
Mpc^ 



= Ai4 X lO-^E{zy 



M200 



IQl^/z-lMcT) 



(9) 



where the pedex indicates that we are considering quantities com- 
puted inside the volume enclosed by /?2()0- By fitting the data com- 
puted from their simulated clusters sample. , Nagai (2006) obtained 
Ai4 = 2.56 and a = 1.70 (we consider the results of their CSF sim- 
ulation, which accounts for a variety of physical processes of the 
baryonic component): we used this relation assuming, as we did for 
the X-ray modeling, a self-similar redshift evolution, correspond- 
ing for the tSZ effect to s = 2/3. 

We show the results of our logA'-logK"" in Fig.|4] Most of the 
conclusions derived from the analysis of the X-ray results apply 
also for the tSZ effect: the uncertainties in the estimate of Og make 
it very challenging to discriminate between different non-Gaussian 
models. However it is worth to notice that, since the tSZ effect is not 
affected by redshift dimming as the X-ray flux, at higher redshifts 
(z > 1) the number of detections is still significant. For example, as- 
suming the flux limit expected for the SPT surv ey (i.e. 5 mjy at 150 
GHz, corresponding to F ~ 3.4 x 10^^^, see iMaiumdar & Mohd 
l2003h . we predict the possibility of detecting ~ 0.3 high-z objects 
per square degree: given the expected area of 4000 deg^, this leads 
to more than one thousand objects, enough to potentially discrimi- 
nate the ~ 20 per cent difference in the cluster counts predicted for 
the /nl= ±100 models. Note that these deviations are comparable 
to an uncertainty of about ±0.02 in the primordial power spectrum 
normalization Og . The expected number counts from the SPT sur- 
vey in different redshift ranges are reported in Table|2] 



4.1 Integrated properties of the tSZ signal 

Apart from cluster counts, other global quantities like the aver- 
age Compton y-parameter and the tSZ power spectrum can be af- 
fected by the presence of a non-Gaussian signature in the primor- 
dial power spectrum of perturbations. In order to study these ob- 
servables, it is necessary to create and analyse mock maps of the y- 
parameter. Since galaxy clusters are extended sources, we need to 
make further assumptions about the density and temperature pro- 
files of the haloes to model the distribution of the signal on the 
sky. Since the tSZ signal receives a significant contribution also 
from the external regions of galaxy clusters, the modelisation must 
take into account the steepening of the slo pe of these profiles in 
the regions around the virial radius (see, e.g.. lVikhlinin et alj2006l : 
iRoncarelli etai]|006bl). In this context, it is imp ortant to note that 
a classic B-model Jcavaliere & Fusco-Femianol Ll978) would fail 
simply because it does not converge for large distances from the 
centre, even when assuming a decrease of the temperature with a 
polytropic index, as adopted in lAmeglio et al.l (I2OO6). Therefore, 
we start from the suggestion of ^ Ameglio et alj (12006) and slightly 
modify their results by introducing a rolling-(3 polytropic profile for 
the tSZ signal, namely 



j(e)=yo 



i± 



-3Perf(7-l)/2 



(10) 



where 9 is the angular separation from the cluster centre, 9c is the 
angular size of the core radius (assumed to be = 0. 1^200)1 7 = 
1.18 is the polytropic index of the gas. The effective slope pgff is 
given by 

j;±Pint/Pext' 



P, 



eff '- 



-Pe: 



1 



(11) 



where x = 9/9c and (3ext and Pim are the external and internal slopes 
for the density profile, respectively. In this way we can tune these 
two para meters in such a way that this expression converges to the 
profile of lAmeglio et al. I ( l2006h (and to a p-model) in the inner part, 
while in the outer part it steepens to allow the value of the integrated 
y'"' to converge to a finite value (ypext > 1)- We choose Pint = 2/3 
and Pext = 1-3, where the latter value is taken in agreement with 
the analyses of Roncarelli et a l. ( 006b) on the density profiles in the 
outskirts of simulated galaxy clusters. With this choice of parame- 
te rs, the external par t of our y{Q) profile also agrees with the results 
of iHaugboelle et all {2007). A visual comparison of these differ- 
ent profiles for Compton y-parameter is shown in the left panel 
of Fig. |5] while the three profiles converge to the same value at 
the centre, the roUing-P polytropic profile is significantly lower al- 
ready at ^ 39c (~ 0.3i?20())- When considering the distribution of 
the signal (right panel), it is easy to deduce that only with the pro- 
file adopted in this work )""' converges to a finite value for high 
values of 9. It is worth to notice that even adopting this profile, the 
total signal up to ^^200 is only the 93 per cent of the total signal: 
this indicates that using any other shallower profile would lead to 
non-negligible biases. 

Then we used this profile as a weight to distribute the total tSZ 
signal of each halo into the pixels of our maps. Since the adopted 
profile has non-zero contribution even at large scales, we set up the 
integration limit to the 99 per cent of the total value: this means 
integrating up to ~ l.S rjoo. We follow the procedure described in 
IRoncarelli et al.l ( l006ah to smooth the signal of the SPH particles 
and obtain the tSZ map corresponding to each simulated light-cone 
realization. At the end our analysis will be based on a total of 60 
maps from the first simulation set and 700 from the second one. As 



8 M. Roncarelli et al. 




-13 



-12 



13 



-12 




-13 



-12 



-10 



-13 



-12 



Log(Y) 



Log(Y) 



-10 



-9 



Figure 4. As in Fig.[T] but for the logW-logF. The dashed vertical line indicates the SPT flux-limit, i.e. 5 mJy at 150 GHz, con'esponding loY 3.4 x 10 



et al, 2006 {(l = 2/i, 7=1.18) 
•■ (ft.,-2/3, (!..,- 1.3, 7-1.18) 




Figure 5. Lef t panel : comparison between the Compton y-parameter profiles for an isothermal p-model (dashed line), the polytropic p-model proposed by 
lAmeglio et alj |200g) (dot-dashed line) and the rolling-p polytropic profile (solid line) used in this work. Right panel: distribution of the signal as a function 
of the distance from the centre for the same three profiles. 



an example, in Fig.|6]we show the Compton v-parameter maps cor- 
responding to the same light-cone realization, as obtained from the 
/nl= 0, ±1000 simulations of the second set. With these extremely 
large amounts of primordial non-Gaussianity, it is possible to rec- 
ognize the expected behavior: an evidence of more (less) clustered 
signal for positive (negative) /nl, compared to the Gaussian case. 
Notice that the average value of the y-parameter computed over all 
the maps of the Gaussian simulation with the WMAP-1 cosmol- 
ogy is < y >= 9.71 X 10^^, while considering a WMAP-3 cosmol- 
ogy (first set of simulations) this value drops to < y >= 5.22 x 



10 ^, in good agreement with the expected scaling, y oc 



.3,5 



(see, e.g.. lKomatsu & Seliakl2002l : lDiego & Mai umdar"200 4l). The 
last fi gure is also consistent with the results of Roncarel li et al.l 
who, analysing a high-resolution hydrodynamical simula- 



tion JBorgani et al. I l2004l) based on the same cosmological model, 
obtained a value of<y>=1.19x 10^* with about half of the sig- 
nal originated from the WHIM (not considered in this work). No- 
tice that the mean values derived for the models with the WMAP- 
3 cosmology but with /nl= 100 and /nl= —100 are < y > = 



X-ray and SZ imprints of primordial non-Gaussianities 9 




Figure 6. Examples of maps for the tSZ signal (expressed in terms of the Compton y-parameter) for three different levels of primordial non-Gaussianity, 
/ml = —1000,0, 1000 (left, center and right panels, respectively). These maps have been obtained from the second simulation set (WMAP-1 cosmology): they 
are 5.71° on a side with a pixel size of (20 arcsec)^. Notice that the three maps refer to the same light-cone realization. 



5.37 X 10 ^ and < y >= 5.08 x 10 ^, respectively, thus compa- 
rable within few percent to the Gaussian case. 



4.2 Angular power spectrum 

To characterize the statistical properties of the tSZ effect, it is im- 
portant to study its power spectrum at different multipoles, in par- 
ticular focusing on the angular scales at which the tSZ signal is 
expected to dominate the primary CMB anisotropics (£ > 2000). 
For the complete set of maps generated as described in the previ- 
ous Section, we computed the tSZ power spectra, using a method 
based on Fast Fourier Transforms, adopting the approximation of 
flat sky (given the reduced extension of the maps) and assuming the 
RJ frequency limit. Finally, the corresponding averages are shown 
in Fig.[7]for the different models. 

Again, when considering the first set of simulations, the dif- 
ferences between models with various values of /ml are very low 
when compared to the variations related to a change of the cos- 
mological scenario: about 10 per cent in both senses for the /nl= 
±100 models, compared to a factor of about 3 when changing 
and i^n,. This is expected as, given the dependence Q Og, the 
difference in the as choice accounts alone for a factor 2.3, with 
the remaining difference associated to the fim parameter. For these 
reasons, the perspective of constraining /nl seems to be quite de- 
manding without an independent derivation of the main cosmolog- 
ical parameters. 

The presence of a possible non-Gaussianit y in post- 
inflat ionary perturbations has been also claimed by ISadeh et aU 
( l2007i) as a possible explanation of the anomalous values of the 
tS Z power spectrum obt ained by the BIMA experiment. In particu- 
lar, |D^son|etal] j2002l) measured an high value of AT = 16.6^5 ^ 
/jK at £ = 5237 (see, however, the smaller estimate obtained by 
ISharp et al.ll200l with the SZA experiment). According to our re- 
sults, explaining the BIMA results with primordial non-Gaussian 
fluctuations alone (and keeping ag = 0.8) would lead to values of 
/nl unrealistically high (/nl^ 100). On the contrary, if slightly 
higher values of ag are considered, the tSZ power spectrum woul d 
agree within Ict with the result obtained bv lOawson et alj j2002h . 
In this context, a positive value of /ml within current upper lim- 
its, could also contribute to boost the tSZ signal in order to explain 
these measurements. 
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Figure 7. The power spectrum of the tSZ signal in the RJ limit as a func- 
tion of the multipole (. for the two simulation sets. The solid lines refer to the 
first simulation set (average of the 20 light-cone realizations) for the three 
different levels of primordial non-Gaussianity (/nl=0,±100), while the dot- 
dashed lines refer to the second simulation set (average of 100 light-cones, 
with WMAP-1 cosmology) for /nl=0,±1000. The dotted line represents 
the primary CMB signal calculated using CMBFAST (Seljak & Zaldarriag^ 
[l^996) adopting the WMAP-3 cosmology. The diamonds with errorbars (lo) 
represent the measurements of iDawson et al.. QOO^ with the BIMA experi- 
ment. 



5 DIFFERENTIAL REDSHIFT COUNTS 

From the results presented above it is clear that a significant de- 
tection of the signatures of a possible primordial non-Gaussianity 
based on global properties of galaxy clusters (e.g. number counts, 
mass functions, etc.) appears very difficult and well beyond the pos- 
sibility of current and planned surveys. The main reason is not only 
the degeneracy of the results with other uncertain cosmological pa- 
rameters (mostly Og), but also the fact that these observables are 
dominated by low-redshift haloes (z < 0.5) where the LSS proper- 
ties of the different non-Gaussian models show smaller differences. 

For these reasons the most reasonable strategy to break the 
degeneracy between /nl and the other parameters can be a study of 
the evolution with redshift of the cluster counts. In Fig. [8] we show 
for the different models under analysis the redshift distribution of 



10 M. Roncarelli et al. 



the objects that will be detected by the eROSlTA and SPT cluster 
surveys (left and right panels) in the X-ray and tSZ, respectively. 

Looking in more detail at the eROSITA results, the cluster 
abundances at z > 0.5, where the satellite is expected to detect 
about 15,000 objects (see Tab. [TJ, show a relative difference of 
about 10 per cent. Therefore, in principle, the possibility of esti- 
mating, thanks to a dedicated follow-up campaign, a high number 
of redshifts for the objects detected by eROSITA at z > 0.5, would 
allow to obtain their redshift distribution, increasing the chance of 
constraining /nl> once the value of is derived using the abun- 
dance of low-redshift objects. 

A similar conclusion can be drawn by looking at the redshift 
distribution of the SPT clusters (right panel). The possibility of de- 
tecting very high-redshift (z > 1) clusters with an SZ survey is cer- 
tainly promising, since in this redshift range the relative differences 
grows to ~ 20 per cent. However, the possibility of obtaining a sig- 
nificant amount of redshift estimates for these objects (which are 
~ 1000, see Section|4j is of course much lower, due to their lower 
signal. 



6 CONCLUSIONS 

In this work we used a set of cosmological W-body simulations to 
investigate the impact of primordial non-Gaussianity (parametrised 
in terms of /nl) on the LSS. From their outputs we constructed 
halo catalogues at different redshifts and, making use of suitable 
scaling relations between masses and observables, we computed 
their expected X-ray emission and SZ signal. In particular we in- 
vestigated the possibility of constraining /nl with future projects, 
namely eROSITA and SPT cluster surveys. Moreover we discussed 
the degeneracy with other uncertain cosmological parameters, like 
Dm and Og. Our main results can be summarised as follows. 

(i) As predicted b y analytical models (see, e.g., lMatarrese et al] 
l200d :l Lo Verde et al . 2008), the main effects in cluster counts due 
to the presence of some level of primordial non-Gaussianity are 
for high masses and redshifts. In particular, for haloes with M200 > 
IO^'^/i^'Mq at z > 1, the differences with respect to the correspond- 
ing Gaussian models are about 20 per cent. 

(ii) When the power spectrum normalization suggested by 
WMAP-3 is adopted, models with a moderate level of non- 
Gaussianity (/nl=±100) well reproduce the observed cluster 
counts derived from ROSAT cluster survey. However, the depen- 
dence on /nl is very weak, when compared to the one on Og, which 
must be independently estimated to fully exploit cluster counts as a 
probe of primordial non-Gaussianity. 

(iii) We predict the expected number and redshift distribution 
of the galaxy clusters that will be detected in two future cluster 
surveys: eROSITA (X-ray) and SPT (SZ). The effects due to a mod- 
erate primordial non-Gaussianity are in general of few per cent, 
reaching about 20 per cent only at high z. In general, the fact that it 
is easier to detect high-z objects with SZ observations, because of 
the absence of cosmological dimming, makes SPT a more promis- 
ing probe for obtaining constraints on /nl- However, once again, 
the results show a strong degeneracy between /nl and other cos- 
mological parameters. Similar conclusions can be also drawn when 
analysing the power spectrum of the tSZ signal produced by galaxy 
clusters. 

(iv) On the whole, the best strategy to detect the signatures of 
primordial non-Gaussianity in the LSS is to perform deep cluster 
surveys, together with suitable optical follow-ups for the determi- 
nation of their redshifts. With this kind of observational dataset, it 



would be possible to constrain ag using low-redshift objects and 
analyse the dN/dz in the range 0.5 z < 1 to constrain the value 
of /nl- If with this method future surveys will allow to reduce the 
uncertainties on Og to about 0.01, this would make possible to de- 
tect moderatejion-Gaus^ order of /nl= ±100 (see 
also lSefusatti et al.ll2006l : rOgurill2009l) . 

In conclusion, the results of this paper confirm the power of 
statistical tests based on galaxy clusters as a probe for primor- 
dial non-Gaussianity. In particular the detection of objects in the 
high-mass tail at sufficiently large redshift, as possible in future SZ 
wide surveys like SPT, will be certainly useful to improve the con- 
straints on /nl coming from alternative methods, like CMB, ISW 
and galaxy biasing. 
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